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We consider the theory obtained by adding to the usual string frame dilaton gravity action specially 
constructed higher derivative terms motivated by the limited curvature construction of |Q] and 
determine the spatially homogeneous and isotropic solutions to the resulting equations of motion. All 
solutions of the resulting theory of gravity with these symmetries are nonsingular and all curvature 
invariants are bounded. For initial conditions inspired by the pre-big-bang scenario solutions exist 
which correspond to a spatially flat Universe starting in a dilaton-dominated superinflationary phase 
with H > and having a smooth transition to an expanding Friedmann Universe with H < 0. Hence, 
the graceful exit problem of pre-big-bang cosmology is solved in a natural way. 
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I. INTRODUCTION 

One of the paramount problems manifest in modern 
cosmological theories is that of the initial singularity. 
Gasperini and Veneziano have proposed the pre-big-bang 
(PBB) model based on the assumption that the Uni- 
verse evolved out of a ten-dimensional Minkowski space- 
time which is an exact vacuum of string theory, in hopes 
of resolving the singularityjproblem. Recently, Branden- 
berger, Easther and Maia O managed to construct non- 
singular cosmological solutions based on initial conditions 
resembling those of the pre-big-bang scenario. They be- 
gan with dilaton gravity and the low energy effective ac- 
tion of string theory in the Einstein frame, adding to 
it specially constructed higher derivative terms. As the 
string frame is, the fundamental frame with respect to 
string theory u , it is important to reconsider the analy- 
sis of Q , performing the calculations directly in the string 
frame. The Einstein frame metric g^n, is related to the 
s^ing frame metric g^^i, via a conformal transformation 



(1) 



where if is the dilaton. 

The (PBB) model is constructed from the low energy 
effective action of string theory. In the string frame, after 
conjpactification to four dimensions, the action is given 
by i 



S 



where = 8ttG = 87rm, 



{R+{V^f + ...}, (2) 



■pi 



The PBB model gives us hope of finding a nonsingu- 
lar cosmology in which in the Einstein frame the Uni- 
verse starts out in a cold dilaton-dominated contracting 
phase, goes through a bounce and then emerges as an 



expanding FRW Universej^] In the string frame this cor- 
responds to a Universe starting in a superinflationary 
phase with H > and emerging as a usual FRW Universe 
with H <0. 

As is well known, it is impossible to smoothly con- 
nect the contracting and expanding branch^pf the PBB 
model using only the tree-level action O^EiJ. These 
branches are separated by a future and past singularity 
of the branches respectively. However, one-loop effects 
in superstring-£osmology are capable of regulating the 
singularities 1^ and smoothly connecting a contracting 
phase to an expanding phase in the presence of spatial 
curvature. Rcfs p^-|2^ describe other attempts to regu- 
late the singularities of pre-big-bang cosmology. Most of 
these approaches, however, have the drawback of being 
perturbative in nature, and that their consequences are 
only felt in the region of the theory which can only be 
described by the full non-perturbative physics. We wish 
to discuss a mechanism to solve the graceful exit problem 
which involves non-perturbative physics. In the absence 
of a non-perturbative definition of string theory we will 
discuss a toy model which only contains gravity and the 
dilaton. 

It is natural to assume that any effective theory of grav- 
ity obtained from string theory, quantum gravity, or by 
integrating out matter fields, will contain higher deriva- 
tive terms in the action. The question we wish to address 
is if by adding specially chosen higher derivative terms 
to the Ricci scalar of the Einstein action we are able to 
cancel the divergences of classical gravity present at ex- 
tremely large curvatures. 

Hence, our general approach will be to construct a 
higher order theory of gravity admitting cosmological so- 
lutions which are everywhere nonsingular and compatible 
with the PBB scenario. 



*See http://www.to.infn.it/'~gasperin/ for an updated col- 
lection of papers on pre-big-bang cosmology. 
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One such construction is carried out in Refs |l|,p9[ . 
The resulting action includes a particular combination of 
quadratic invariants of the Riemann tensor, similar to a 
Gauss-Bonnett term added to the usual Einstein- Hilbert 
action for gravity. The inclusion of this invariant causes 
all solutions of the equations of motion to approach de 
Sitter space-time at high curvature, which is naturally 
nonsingular. Although in general most higher-derivative 
theories of gravity have much worse singularity properties 
than Einstein gravity, this particular construction stems 
from the "limiting curvature hypothesis" (which we will 
discuss more below) and is entirely nonsingular. 

The most elementary way to formulate a nonsingular 
gravitational theory is to add an invariant I2 to the ac- 
tion, with the property that /2 = is true if and only if 
the spacetime is de Sitter. By coupling I2 into the grav- 
itational action via a Lagrange multiplier field ip with a 
potential chosen to ensure that /2 ^ at large curvatures 
we can force all solutions to approach de Sitter space in 
this regime. For homogeneous and isotropic spacetimes, 
a simplepChoice for I2 which satisfies the aforementioned 
criteria is 



(3) 



Recall that the physically measurable curvature invari- 
ants such as R and R^j^^R^'^ are constructed from the 
string frame metric, g^v The simplest desired action 
including I2 is 



S{g^i,,,i') 



\R + i;l2-V{4>)) (4) 



where V{ip) is a function chosen such that the action has 
the correct Einsteinian low curvature limit for —^ 0, 
whereas ioi ^ 00 the constraint equation forces I2 
0. After integrating out -0, we obtain a higher derivative 
gravity theory which is nonpolynomial in the curvature. 

Note that introducing the Lagrange multiplier field 
is a convenient way of writing a non-polynomial higher 
derivative gravity action in a polynomial form. ip has 
no independent physical meaning. By itself, it neither 
represents a new matter field, nor does it stand for a 
particular symmetry. The situation is analogous to the 
construction of the action Sr = J dt\/l — for point 
particle motion in special relativity from the correspond- 
ing action S^n — J dt^x^ for point particle motion in 
Newtonian mechanics. If the Lagrange multiplier con- 
struction analogous to the one employed in this paper is 
used to construct the new theory, with a general poten- 
tial satisfying the asymptotic conditions required, then 
the resulting theory will have bounded velocity but will 
not be the theory for point particle motion in the Special 
Relativity, i.e. in the theory with the new symmetry. Sji 
can only be obtained by means of a very special choice 
of the coefficients of the potential for ip. In analogy, we 
cannot expect to be able to extract any new stringy sym- 
metry from the ad-hoc construction of this work. 



In this paper, we will show that the addition of the 
same higher derivative terms as used above in (^) to the 
action of PBB cosmology results in the elimination of sin- 
gularities and connects the (in the Einstein frame) dila- 
ton driven contracting phase with the expanding phase 
via a bounce. 

From the point of view of string theory our construc- 
tion appears, unfortunately, rather artificial in that we 
are picking out a particular narrow class of (nonpolyno- 
mial) higher derivative gravity actions from the infinite 
dimensional space of such actions. Our work should be 
viewed as an existence proof (in the string frame) of a 
higher derivative dilaton gravity theory with nonsingular 
cosmological solutions rather than as a realistic string- 
motivated model. 

In particular, the invariant I2 which plays a crucial role 
in our construction is not one of the curvature invariants 
which appears in the next to leading order in derivatives 
of the low energy effective action of string theory. The 
role of the invariant I2 and of the associated Lagrange 
multiplier field is to implement the limiting curvature 
hypothesis. There are indications that this hypothesis is 
realized in nonperturbative string theory (see e.g. pof), 
as a consequence of t-duality. This duality is broken when 
considering the low energy effective action of string the- 
ory, and it is therefore not to be expected that the cur- 
vature invariants which arise in this low energy effective 
action can be used in place of 12- The invariant I2 plays 
the role of non-perturbative information about the the- 
ory. Obviously, since the non-perturbative definition of 
string theory is still unknown, we cannot hope that I2 
yields more than a toy model for what string theory will 
eventually tell us. Since the limiting curvature hypothe- 
sis plays an important role in both our construction and 
in string theory, we can hope that our model may even- 
tually turn out to contain some of the relevant physics. 

We should also warn the reader that our work has poth- 
ing to say concerning the initial condition problem of 
the pre-big-bang scenario. 



II. ACTION AND EQUATIONS OF MOTION 

We begin with the string frame action for dilaton grav- 
ity, adding to it the higher derivative term given by I2 , in 
analogy to what was done in the absence of the dilaton 
in Refs 0j2^ and in the Einstein frame in Ref 

^ " 2^ /'^^^ N/^e-"^ {R+ (V(^)2 + c^e-'^^h - Vi^)}. 

(5) 

For the moment, we allow a general coupling between I2 
and the dilaton. Minimal coupling corresponds to setting 
the constant 7 equal to zero. The constant c rescales the 
Lagrange multiplier field tp, and will be chosen to simplify 
the equations of motion. 
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Consider a homogeneous and isotropic metric of FRW 
type 



1 



{t)dt^ - a{tf{ j-^dr^ + r^dfl^) , (6) 

1 



where dfl^ is the metric on S'^ and n{t) is an arbitrary 
lapse function which we wiU set to unity in the EOM 
below. The equations of motion resulting from varying 
with respect to (p, ip and n become 



H 



-If dV 



cVV2 dtp 

sT"^\/T2 (sH^iP --^tP + Hi^jj + (7 - l)iptp) 



(7) 



respectively, where dots denote differentiation with re- 
spect to time, t. 

For a spatially flat, bouncing Universe, we set the cur- 
vature constant fc = 0. We will simplify things further by 
considering minimal coupling of ip to I2, i.e. 7 = 0. To 
eliminate useless constant coefficients in the equations 
of motion, it is convenient to choose c = 17-^12. The 
resulting equations of motion become 



V' = 

H = 
X = 



6H +{x-6)^J + -^(x' 



1 



(8) 



with X — V' ^-i^d the prime (/) denoting differentiation 
with respect to ip. 

Now let us focus on the construction of the potential 
V{ip). Once again, our arguments are the same as in Ref 

t. When the curvature is small, the terms in the action 
with Lagrange multiplier field tp dependence must be 
negligible compared to the usual terms of dilaton gravity. 
This is made manifest by the condition 



V{tP) ~ \tp\ 



(9) 



since the region ofj-small \tp\ will correspond to the low 
curvature regime. O In order to implement the limit- 
ing curvature hypothesis, the invariant I2 must tend to 
zero, and hence the metric ij^jy will approach a de Sit- 
ter metric at large curvatures, i.e. for 1-01 — > 00. From 
the variational equation with respect to tp, we find the 
requirement 



const \tp\ 



(10) 



As in ^ we add a third condition in order to ensure 
that there is a bouncing solution for k = 0. The equa- 
tions must allow a configuration with H — and tp 0. 



From the equation of motion for tp in (and considering 
the region of small x) it follows that V{tp) must become 
negative, assuming that it is positive for small \tp\. Let 
tpf, denote the value of tp at the bounce. This will occur 
when 



(11) 



as can be seen from the equation of motion for tp setting 
H = 0. 

A simple potential which satisfies the conditions (|^), 
(0) and goes negative beyond some critical value of tjj is 



V{tP)^ 



(12) 



Note that the specific values of the coefficients of the 
terms in V{tp) are not important as long as the three 
criteria discussed above are satisfied. 



III. PHASE DIAGRAM OF SOLUTIONS 



The condition (11) on the potential V{tp) ensures that 
our model is nonsingular and geodesically complete for 
large values \tp\ but says nothing about geodesic com- 
pleteness of solutions which always remain within the 
small \tp\ regime. Hence it is necessary to study the small 
I ■01 region in order to prove that our model is everywhere 
nonsingular. 

In this section we discuss the projection of the three 
dimensional phase space {tp, H, x) of our model into the 
two dimensional phase plane {tp, H) for small values of x- 
By projecting out the x axis we lose no relevant informa- 
tion about the system since the curvature invariants are 
not affected by x- Both analytical and numerical meth- 
ods are used to study the trajectories of solutions of (H) 
in the phase plane and thus explicitly demonstrate the 
absence of singularities. Furthermore, and unlike in the 
model of Refs. ||,|2^, we shall show that our theory ad- 
mits spatially flat bouncing solutions as in the Einstein 
frame model of 

We need to demonstrate that all solutions tend to a 
known nonsingular and geodesically complete space-time 
at large positive and negative times. With the symme- 
tries of our problem, these asymptotic space-times are 
either Minkowski or de Sitter. We need to show that all 
phase space trajectories either tend to the origin of phase 
space or else to \tp\ 00 for finite values of H. We will 
first analyze the phase space trajectories for small values 
of |x|. In Section 4 we then discuss the case of large |x|- 

There are many interesting points and curves on the 
phase plane {tp, H) which demand special attention. 
First, the point {tp,H,x) = (0,0,0) corresponds to 
Minkowski space-time. The potential V{tp) vanishes at 
this point, but it also vanishes at the values 



0b = ±4 . 



(13) 
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From the first equation in (||) it follows that the phase 
plane points (V'bjO, 0) correspond to bouncing points of 
cosmological trajectories. The general equation for a pos- 
sible bounce is 



(14) 



To prove that this is in fact a bounce, we need to show 
that contracting solutions are attracted to it. This will 
be done below. 

H vanishes when the derivative of V does, i.e. at the 
values 



^rf = ±2 . 



(15) 



Trajectories which cross the phase space planes {tpd, H, x) 
have ij = at the points where they cross. 

To show that the curves through {ip, H, x) = (4 — 
/(x),0,x), with /(x) chosen such that (x^ - V) ^ 0, 
represent bounce solutions, we expand the ip equation of 
motion near H = 0, which gives 



(16) 



First note that contracting solutions with ^ 1 and 
2 < -0 < 4 — /(x) have ip > {) and approach the bounce 
(4— fix) J Oj x) in finite time since H is positive and does 
not go to zero. If the trajectories reach ip > i — /(x) 
before hitting the ip axis, then -0 changes sign. This 
shows that the bounce is an attractor on the contracting 
branch. By taking the time derivative of the ip equation 
in (^ and using the other two equations to eliminate x 
and H and after expanding near the bounce it can be 
verified that on the expanding side of the branch ip < 
and that therefore the trajectories will eventually turn 
around (i.e. tp < 0) (see Figure 1). 

By expanding the equations about the origin we find 
that 



d^P 1 



6H^ + x' 



which defines the critical surface 



0^ + GiPH - 6H^ - X 



6H) 







(17) 



(18) 



Note that lines on the critical surface are not phase space 
trajectories of the equations of motion. Such trajectories 
would point in the vertical direction at the critical surface 
since dip/dH = 0. Note also that as x^ increases, the 
critical surface approaches the ■0 axis. 

For the contracting solutions we find that the equation 
of motion for H becomes 



H 



-2^p. 



(19) 



Hence contracting solutions with -0 > which are above 
the critical surface have H < and from Eq. ([l^) -0 > 0. 
These trajectories are pulled toward the tp — 2 curve 
where H changes sign. Trajectories which start out in 



0.2 



. 1 



-0 . 1 



-0.2 




FIG. 1. Projection of a set of trajectories onto the tp - 
H plane of phase space for small values of and X- The 
Hubble parameter, H, is plotted on the vertical axis, while 
tp is plotted on the horizontal. Shown are a set of solutions 
which oscillate about Minkowski space-time H = 0, tp = 
as well as solutions through the bounce points H — and 
tp ~ ±4. 



this region, do not cross the critical surface and reach 
the = 2 surface at a value \H\ <C 1 are good can- 
didates for spatially flat bouncing Universes (see Fig- 
ure 1). The analysis of trajectories through {tp,H,x) — 
(—4 + /(x))0,x) is similar. Trajectories near the origin 
below the critical surface correspond to solution which 
oscillate about Minkowski space-time (see Figure 1). 

There is a separatrix surface in the lower right quad- 
rangle of the phase plane which separates trajectories 
which reach tp — 2 from those which do not and which 
eventually cross the critical line. The solutions for tra- 
jectories above this separatrix line close to the origin can 
be obtained by taking the time derivative of the second 
equation in (@) and inserting the first equation and by 
considering terms which dominate near the origin. For 
small values of x^ we obtain 



Hii = -6HH + , 



(20) 



which is the same equation as the corresponding one in 
except for the first term on the right hand side. However, 
by inserting as an ansatz the solutions found in H 



H{t) = ~ct^ 



(21) 



where c is a constant, it is not hard to show that for small 
values of i, the extra term in our equation gives negligible 
corrections to these trajectories. For sufficiently small 
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values of c, the trajectories lead to bouncing solutions, 
whereas for somewhat larger values of c, the trajectories 
cross the ■0 = 2 plane at \H\ > 1. 

For the latter solutions, the ip equation of motion for 
tp ^ 2 becomes 



^ = _(3F + 6)7^ + lx' 



(22) 



For large values of and not too large values of x^, "0 > 
and H ~i and the trajectory tends to de Sitter space. 
However, as will be discussed in Section 4, x is rapidly 
growing, and eventually (unless the dilaton is stabilized), 
the x^ term in the above equation takes over and leads to 
i/) < 0, in which case these solutions also tend to de Sitter 
space, but at large negative values of tp. The bottom line 
is that there are no singularities in the phase space region 
> 2 and < 0. 

From Figure 2 we see a potentially troublesome fea- 
ture of the phase space near ■0 = 2. There are critical 
lines with ■0 = which converge to ^0 = 2 and H = oo 
from the left, and to ^0 = 2 and H = —oo from the right. 
Along these critical lines, the trajectories appear to head 
off to infinity (for negative H) or emerge from infinity 
(for positive H), suggesting a singularity is present in 
the model. Upon a more detailed analysis of the EOM 
it becomes clear that the solutions in this region are re- 
pelled by the critical line. To see this begin with the 
critical line obtained by setting ^0 = 0. Considering the 
lower right quadrant of the phase plane, to the left of the 
critical line, we have 



6H^ +X^-V 
H{3H + 6-x) 



(23) 



The critical line for small values of approaches the 
surface tp — 2 from the right (note that in the Einstein 
frame model of it approached from the left, the dif- 
ference being due to the extra term in the equation of 
motion in (g)). Differentiating the ip equation of motion 
with respect to time and evaluating the result along the 
critical line we see that 



ip <0 



(24) 



and hence the critical line is a repeller, sending the tra- 
jectories away from the line and toward the asymptotic 
de Sitter region. To the right of the critical line > 
and the trajectories are again pushed away from ^0 = 2. 

Consider the example of a frozen dilaton, x ^ ^■ 
Hence, for \H\ >> |0| we find 



i> ~ {V + 3H^){'4> - 2) 



(25) 



and ^0 < as expected. Solutions to the left of the critical 
line peel away to the de Sitter region while those to the 
right of the critical line approach from de Sitter. 

We also wish to show that \H\ oo if and only if 
■0-^-00. For this we will consider the specific example of 




FIG. 2. Projection of a set of trajectories onto the ip - H 
plane of phase space for larger values of \H\. The axes are 
the same as in Fig. 1. The only regions in phase space where 
trajectories tend to \H\ — oo, i.e. to curvature singularities, 
are near the two critical lines which tend to = 2. As shown 
in the text, these lines repel trajectories, and hence there are 
no singular trajectories even in these regions. 



the small negative ^0, large negative H behavior. In this 
region the change in ip with respect to H is given by 



dip _ 3iJ(0 - 2) 
dH ~ V' 



(26) 



Hence we see that the slope of the solutions increase and 
begin to level off as \H\ increases. 

To conclude this section, we have shown at this point 
that there are no singularities for small values of and 
that all solutions can be continued to arbitrary large 
proper time, which demonstrates geodesic completeness. 
In this following section we turn to the discussion of the 
dynamics for large x^- 



IV. DILATON EVOLUTION AND 
STABILIZATION 

The dilaton is the massless scalar field with gravita- 
tional strength couplings, found in all perturbative string 
theories. It is generally believed that dynamical effects 



will generate a mass 
ken supeijsymmetry 



ii the dilaton in vacua with bro- 
If the dilaton were to remain 
massless & it would affect the values of gauge couplings, 
and produce potentially observable consequences. 

In our model there are singularities in x- Consider the 
EOM for X in the large x^ limit 
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X 2^ 



This has solution 



x{t) 



which blows up when 
x{to) 



(t - to) 



(27) 



(28) 



(29) 



This singularity does not concern us however, since 
it can be avoided by adding a simple potential U{ip) to 
the dilaton equation in order to freeze the dilaton in a 
manner consistent with the above argument. Where this 
potential comes from is the usual problem of any theory 
with a dilaton, a problem to which we have no answer 
either. 

What happens to the projection of the phase space 
trajectories onto the ^ — H plane in the large x regime? 
Here we have 



H 



(30) 



and the H equation of motion remains the same. Imme- 
diately we see that for 0, |-0| ^ cxd which implies 
that '0 — !■ oo and hence H — s- const. Further insight is 
gained by noticing that 



■0 > for H > 
< for H < 0. 



(31) 
(32) 



We know that V vanishes for ijj — Q and = ±2. Thus 



= 0: ^ = ^ 
V' = ±2: = ±2x4 



x_ 

H 



(33) 
(34) 



Therefore, solutions above H — stretch straight (more 
or less) across the ^ — H plane from left to right, and 
those below the H = line stretch from right to left, as 
can be seen in Figure 3. 

In the large x^ limit, the equation for the critical sur- 
face -0 = becomes 



X 

3ff2 



(35) 



However, since for any solution of the equations of motion 
in this regime x'^ grows much faster than 0, trajectories 
do not reach the critical surface. Instead, H ^ 0. This 
completes the proof that there are no curvature singular- 
ities for large x^, and that the solutions asymptotically 
approach either Minkowksi or de Sitter. 

Finally, let us assume that the dilaton has been fixed. 
In this case, x — and the remaining equations of mo- 
tion become almost identical to the ones studied in |^] 
(for X = 0) except for the —6ip term on the right hand 




FIG. 3. The phase portrait for solutions of the equations of 
motion for fixed dilaton. The axes are the same as in Figs. 1 
and 2. Due to the extra terms in the equations of motion, the 
phase portrait is not symmetric under reflections about the tp 
axis as it was in |^ . Clearly visible are the three critical lines. 
As discussed in the text, there are no singular trajectories 
even in the vicinity of these critical lines. 



side of the equation of motion in (|^) . This term, how- 
ever, is subdominant in the large \H\ region, and hence 
the conclusions of concerning the absence of singulari- 
ties carry over to our case. The extra term does, however, 
lead to asymmetries in the phase portrait under reflec- 
tions about the ip axis. There are now three critical lines 
where -0 = 0. The first line occurs for < 0^ < 2 for 
positive values of H, and with _ff — > oo as 0" — > oo. The 
second line occurs for 2 < -0 < oo and for negative val- 
ues of H with H — oo as ^0 — > 2. Finally, for -0 < 
there is a harmless critical line which remains at small 
negative values of H (see Figure 3). As in Section 3 it 
can be shown that the two first critical lines repel the 
trajectories. Hence, the proof of non-singularity of the 
model carries over. 



V. CONCLUSIONS 

We have extended the model of nonsingular dilaton 
cosmology presented in to the string frame. Us- 
ing the limiting curvature construction of [l|,^9[| applied 
to dilaton gravity formulated in the string frame, we 
have obtained spatially flat bouncing cosmological so- 
lutions. The construction consists of adding specially 
chosen higher derivative gravity terms in the form of a 
curvature invariant I2 to the string frame action. The 
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invariant I2 is made up of invariants quadratic in the 
Riemann curvature and has the property that /2 = sin- 
gles out the maximally symmetric de Sitter space-times 
among all homogeneous and isotropic solutions. It is cou- 
pled to dilaton gravity via a Lagrange multiplier field ip. 
The V' field is nondynamical but has a potential V^(V')i 
which was chosen to allow non-singular bouncing solu- 
tions. 

The three dimensional phase space of trajectories was 
studied both analytically and numerically to demonstrate 
that all solutions are nonsingular. Specifically we studied 
a large class of solutions which lead to bouncing cosmolo- 
gies. The dynamics of the bounce are governed by the 
higher derivative gravity terms introduced by the limit- 
ing curvature construction. 

The connection with pre-big-bang cosmology appears 
in a different form than in the Einstein frame. In the 
string frame picture of pre-big-bang cosmology, the Uni- 
verse starts in a superexponentially expanding dilaton- 
dominated phase with H > and H > 0, i.e. in the 
upper left quadrangle of the projected phase space of Fig- 
ures 2 and 3. In the absence of the higher derivate terms, 
the trajectories would diverge to = 00. However, as 
is obvious from Figures 2 and 3, the new terms we have 
added lead to a graceful exit from this phase. The tra- 
jectories cross the 7/1 = axis and go over to trajectories 
which are like the usual expanding FRW solutions with 
H > and H < 0. This happens independent of whether 
the dilaton is frozen at late times or not. 

There is a singularity in the dilaton equation of motion 
but it is assumed that this problem can be solved by the 
introduction of a dilaton potential U{ip). Although we 
do not propose a specific form for this potential, such 
a potential is physically necessary in all theories with a 
dilaton in order to generate a mass for the dilaton. 

The most obvious criticism of this model is that the 
higher order terms in the action are artificially con- 
structed rather than derived from fundamental physics. 
However it is important to recall that our method is well 
motivated. For one, all effective theories of gravity, in- 
cluding those produced by string theory, quantum grav- 
ity, or from quantizing matter fields in curved spacetime 
must contain higher derivative terms in the action. Fur- 
thermore, it is natural to assume that physical invariants 
must be limited in such theories to avoid singularities. 
We have found that in both the Einstein frame and the 
string frame, the method proposed here ensures that all 
physical invariants remain finite, and produces bounc- 
ing nonsingular cosmological solutions. The graceful exit 
problem of pre-big-bang cosmology is solved naturally by 
our construction both in the Einstein and string frames. 
The connection with pre-big-bang cosmology will be ex- 
plored further in a future publication. 
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